Hamiltonian theory of the strongly-coupled limit of the Kondo problem in the 

overscreened case 



Domenico Giuliano ^ and Arturo Tagliacozzo * 
' Dipartimento di Fisica, Universita della Calabria, and 
I.N.F.N., gruppo collegato di Cosenza 
Arcavacata di Rende ( CS), 1-87036, Italy 
* Coherentia - (INFM), Unita di Napoli, Napoli, Italy, and 
Dipartimento di Scienze Fisiche Universita di Napoli "Federico II ", 
Monte S.Angelo - via Cintia, 1-80126 Napoli, Italy 
(February 2, 2008) 



By properly generalizing Nozieres' Fermi liquid theory, we construct an Hamiltonian approach to 
the scattering of conduction electrons off a spin-f/2 impurity in the ovescreneed Kondo regime, 
as T — » 0. We derive the S-matrix at the interacting fixed point, and the corresponding phase 
shifts, together with leading energy corrections to the unitary limit. We apply our results to obtain 
the low-temperature dependence of the 2-channel Kondo conductance, and we relate it to possible 
transport experiments in a Quantum Dot. 
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I. INTRODUCTION 



The Kondo Effect in metals containing magnetic impurities consists of an "anomalous" minimum in the resistivity 
p(T), as the temperature T drops below the "Kondo temperature" Tk- The minimum is due to antifcrromagnctic 
scattering of conduction electrons off the localized magnetic impurities [1]. Quite recently [2,3], the signature of 
Kondo interaction has been found in transport experiments across a Quantum Dot (CD) at Coulomb Blockade (CB). 
At CB, the QD is expected to be insulating, due to the discreteness of its levels [3]. Nevertheless, within a CB valley, 
the linear DC conductance may saturate at low temperatures. As T > Tk, the conductance exhibits the typical 
logarithmic raise [4]. When T drops below Tk, it saturates to its unitary limit, G — 2e 2 /h. Below Tk, a perturbative 
approach in the coupling is not feasible. 

The single impurity Kondo effect is classified according to the spin s of the impurity and to the number of channels 
k of itinerant electrons involved in the scattering. Indeed, electrons scattering off the impurity may be labelled 
by quantum numbers other than the spin (for instance, angular momentum) [5]. In spite of the fact that in the 
perturbative temperature region there are no qualitative differences between one-channel and many-channel effect [6] , 
because of the different nature of the corresponding ground states (GS), deep differences arise in the unitary limit, 
depending on k and s. In particular: 

• If k — 2s, as T is lowered down to 0, the flow of the coupling strength between the impurity magnetic moment 
and the spin of itinerant electrons runs all the way towards an infinite-coupling fixed point. At the fixed point, 
the impurity spin is fully screened, and the localized magnetic moment effectively disappears. The impurity 
rather works as a spinless scattering center in the Fermi sea of the itinerant electrons, and double occupancy 
at the impurity site is forbidden ("Nozieres picture" [7]). This is the most common case: a QD with an odd 
number of electrons behaves as a spin-1/2 impurity, interacting with one channel of conduction electrons from 
the Fermi sea of the contacts. 

• If k < 2s, the system is "underscreened" . A residual magnetic moment survives on the impurity, even at T = 0. 
This moment interacts ferromagnetically with the spins of the itinerant electrons. T = is again an attractive 
infinite-coupling fixed point and, accordingly, the system is described by Fermi Liquid (FL) theory. An example 
of this case occurs in a two-dimensional QD in a magnetic field at a singlet-triplet level crossing [8-10]. 

• If k > 2s and the exchange coupling is the same for both channels, the impurity magnetic moment gets "over- 
screened" . An effective residual magnetic moment survives on the impurity, antifcrromagnetically interacting 
with the spin of itinerant electrons. As a consequence, the infinite-coupling fixed point is now repulsive, as well 
as the noninteracting fixed point. The system flows toward an intermediate, finite-coupling fixed point, where 
FL theory has been predicted to break down [11]. 

The prototype model for this case is the two-channel spin 1/2 Kondo (2CK) model, whose low-temperature 
behavior will be the subject of this paper. Our approach extends Nozieres' picture of the one channel Kondo 
fixed point [7] to the 2CK case. This allows us to highlight the Non Fermi Liquid (NFL) properties within a 
scattering matrix (S— matrix ) approach, which is particularly suitable to studying the conductance across a 
QD. 

Overscreened Kondo effect has been hypothized to be the driving mechanism of low-energy physics of several 
systems, although, depending on the physical system involved, the relevant degree of freedom of the impurity may 
not be its spin, but some orbital angular momentum labeling its energy levels ("orbital Kondo") [12]. 

For instance, overscreened Kondo effect has been predicted to possibly take place in glassy metals [13]. Whether 
it really arises in these systems, or not, is still a debated question [14]. A two channel Kondo behavior has been 
invoked in an experiment by Ralph and Buhrman on clean Cu point contacts [15]. A different route toward the 
realization of two-channel Kondo effect in a controlled way has been recently proposed in vertical QD's at Coulomb 
blockade [16], or in similar mesoscopic devices [17,18]. According to the predictions concerning single impurity Kondo 
models, as T < Tk, the temperature dependent corrections to the conductance should display a crossover from log- like 
dependence on T to a dependence on (T/Tk) 2 for the perfectly screened case, on ^JT/Tk for the overscreened case. 
In fact, the crossover to a T 2 -dependence has been experimentally seen in a vertical QD with few electrons [4]. No 
experimental evidence for 2CK in dots has been produced yet. 

On the theoretical side, due to its connections to many condensed matter problems, the 2CK model has been 
studied with a number of different techniques [19]. Finite- T corrections have been derived by means, for instance, 
of Bethe-ansatz like exact solutions [20]. Ludwig and Affleck applied Conformal Field Theory (CFT) techniques to 
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determine finite-T corrections to the unitary limit, Wilson ratios and several exact results concerning Green's functions 
[11]. Numerical Rcnormalization Group (RG) techniques applied to multichannel Kondo have a long history [21-24]. 
Abelian bosonization [25] and subsequent refermionization [26] has been used as well as Majorana fermions [27], to 
study the strongly coupled states, by removing the unscattered degrees of freedom. Functional integral methods, 
including the slave bosons technique [28] and the Coulomb gas approach, have also been applied to 2CK [29]. 

A number of approaches have been applied to describe nonequilibrium properties such as the nonlinear conductance, 
mostly in connection with the Anderson impurity model [2]. These range from numerical rcnormalization group, to 
non-crossing approximation (NCA) [30,31], perturbative functional integral methods [32], perturbative renormaliza- 
tion group methods in real time [33] . 

In this paper we study the 2CK model close to the fixed point, by applying bosonization and refermionization of the 
quantum particle fields. We consider a spin 1/2 impurity at x = 0, antiferromagnetically coupled to a right (R) and 
a left (L) one dimensional noninteracting Fermi sea, with an extra index a = 1,2, which labels the two conduction 
channels. We employ a scattering approach, that is appropriate to study the unitary limit of the conductance at 
T = 0. In particular, this will allow us to calculate both fixed point conductance, and the leading temperature 
dependent corrections. By removing the degrees of freedom not interacting with the impurity in the unitary limit, we 
move the NFL fixed point to infinite coupling. Accordingly, we apply a perturbative strong coupling expansion. We 
first derive the properties of the GS at the fixed point in the bosonic representation using the lattice version of the 
model. Next, going back to fermions in the continuum limit, we obtain a fermionic representation for the fixed point 
S matrix, and get the unitary limit of the conductance, by using the Landauer formula. Finally, we derive finite T 
corrections to the DC conductance, in bosonic coordinates. In the conclusions we propose a connection between this 
model and a possible realization of the 2CK conductance in a QD [16]. 

Our representation of the S matrix does not suffer of the "unitarity paradox" , since, following Ludwig and Malda- 
cena, we introduce a "spin-flavor" quantum number in the bosonic representation [34]. Indeed, in the unitary limit, 
the spin flavor is the only quantum number that is alleged to change upon scattering off the impurity. However, some 
care has to be used when describing such a dynamics in fermionic coordinates. 

The paper is organized as follows: 

In Section II we introduce the Hamiltonian with linearized bands close to the Fermi points and the bosonic rep- 
resentation of the relevant quantum fields. Because of the redundancy, due to the spin- flavor field, we define two 
different bosonic representations for the same fermionic field, which we refer to as / and II. 

In Section III we construct the fixed point impurity GS by reformulating in bosonic coordinates the regularization 
scheme proposed in Ref. [27]. 

In Section IV we go back to the fermionic representation. Once we have identified the physical states, we implement 
Nozieres' scheme, by deriving a one-body potential for the fixed point fermionic Hamiltonian. This rephrases in 
fermionic coordinates the spin- flavor bosonic field dynamics due to the scattering off the impurity. Next, we use 
the Hamiltonian we derive, to calculate the Green's functions at the fixed point. The one-particle Green's functions 
provide us the S matrix elements in the I,i7 representation. Our approach gives the correct result for the phase shift 
for each fermion field, given by ±7r/4. 

In Section V and VI we employ a Schricffcr- Wolff-like transformation, to derive the corrections to the fixed point 
Hamiltonian. In particular, in Section V we show that the first correction, although irrelevant for what concerns the 
fixed point dynamics, selects the appropriate physical states at any point. The unitary limit for the conductance 
follows immediately, provided the degrees of freedom are properly counted. In Section VI we derive the first irrelevant 
operator, giving an energy dependent correction to the phase shifts in the S— matrix and, consequently, a T-dependcnt 
correction to the conductance. 

In Section VII we summarize our conclusions and relate our results to the existing experimental quest for 2CK in 
QD's. 

Mathematical details of the derivation are reported in appendices A,B,C, and D. 

II. THE TWO-CHANNEL KONDO HAMILTONIAN: LOW-ENERGY FERMION MODES AND 

BOSONIZATION 

In this Section we introduce the model Hamiltonian for lead electrons in the overscreened Kondo effect. Since, in 
the following, we will need the lattice version of the Hamiltonian, we start by introducing the lattice version of the 
theory. We will propose the fermionic and the bosonic version of the model Hamiltonian. 

On a system of size L, the lattice kinetic energy in Fourier space is given by: 
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H T = ^ tcos (kea)}cl a {k e )c aa (ki) (1) 

k?;aa 

(the c's are fermionic operator in momentum space, kg = 2n£/L ;£ = ... N — 1, [i is the chemical potential, a is the 
lattice step, a is the channel index). 

In the long wavelegth limit, expanding about the two Fermi points, one gets the effective Hamiltonian: 

H T = -iv f f da;5Z{^ i<ra (ic)^0fl;aa(aj)-0i,; <7a (ic)^0i;<7a(a:)} , (2) 

aa ' 

where Vf = 4tasm(akp)/2ir, and: 

(t>L/R,* a (x) = J dpe ipx (f> L/RtrTa (p) . (3) 

The isotropic lattice Kondo interaction Hamiltonian, is given by: 

H 2 K CK = JS d -[e?i(0)+<?2(0)] , (4) 

where a a (x) — \ a c ta( x )^aa' c a'a( x ) an< ^ is the spin 1/2 impurity located at x = 0. 

By using the linear combinations: 

4>e;aa{ x ) = [4>R- aa (x) + (j>L;aa{-x)] 

<Po;aa{ x ) = -^=[(j)R-aa( x ) - <f>L;<ra(-x)] , (5) 

H]f K takes the form: 



H 2 K CK = JS d • [a e;1 (0) + a e;2 (0)] , (6) 



where: 



that is, only the "e" -fields enter the Kondo interaction Hamiltonian. As a consequence, we may study the Kondo 
dynamics by taking into account only the chiral fields </> e:(JQ . From now on, we will drop the suffix e from the various 
field operators. 

In order to properly deal with the interacting fields, we bosonize (f> aa [11]. Since we have four independent fermionic 
fields, we need the same number of independent bosonic fields, ^ aa - Therefore, following [34], we define: 

4>aa(x)=Va a :e- i9 "-^: , (7) 

where r\ aa are real Klein factors, obeying the anticommutator algebra {r) aa , r\ r {i\ — bar&afi- 

It is possible to introduce the densities of physical quantities starting from the linear combinations [34,26]: 

=^* Q(T (x) ; * sp (a;) = ^ <r* a<r (x) ; * fl (x) = ^ a^ aa {x) (8) 

aa aa aa 

(In Eq.(8) and in the following, whenever we use a and a as coefficients we mean +, —1, when a —1, J,, and +, — 1, 
when a = 1, 2). 

The densities of charge, spin and flavor, p c h, p sp , PR, are given by: 

Pch/ sp /fl(a;) = 2^^* ch / s P/ fl ( x ) ' 

Out of the fields $ M , a fourth bosonic field, the "spin- flavor" field, independent of the first three ones, may be 
constructed, given by: 
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* rf (a;) = ^aa^ aiJ {x) . (9) 
The free dynamics of the bosonic fields *'s is given by the bosonized version of Eq.(2), that is: 

, ■ / dx y \ ) 

J X=ch,sp,fl,sf V 

Notice that the spin-flavor quantum number appears to be "redundant" , as the state of the lead electrons is fully 
determined by charge, spin and flavor. Indeed, it is possible to realize two "inequivalent" representations of the fields 
4> aa in terms of the four fields (A =ch,sp,fl,sf). The former representation, which we will refer to as is given 
by: 

^ a ( X ) = 7? ffQ : e -i[*ch(x)+ CT * sp (x)+a* fl (x)+aa* 8f (x)] . _ (n) 

The latter representation, instead, is defined by: 

<j>™ (x) — : e ~5[* ch ( a: )+ <7 * s p( x )+ a * fl ( x )- Q:CT * 8f ( x )] : . (12) 

The "Klein-like" factors £ CTQ are determined by the requirement that the fields in the two representations anticom- 
mute with each other. Such a requirement is achieved upon defining: 

where the operators N aa are given by: 

N aa = J dx [pch(x) + crp sp (x) + apti{x) - aap si (x)\ 

(notice the unusual definition of the £'s, which involves non real fermionic factors. Nevertheless, both / and II- 
representations provide perfectly legitimate fermionic fields). 

In Appendix A we prove that fields within the same representations obey the usual anticommutation relations, 
while fields from different representations anticommute with each other, that is: 

Ka(*)> <t>%(v)} = 6 ab S aT 6 a p6(x - y) (13) 

(a, b = 1,11). 

We now start the derivation of the effective Hamiltonian in the unitary limit. 



III. FIXED POINT IMPURITY GROUND STATE IN BOSONIC COORDINATES 

To derive the effective theory for the spin- 1/2 overscreened Kondo system in the unitary limit, we use the regu- 
larization scheme introduced in Ref. [27]. Such an approach allows for moving the intermediate coupling fixed point 
towards infinite coupling. In particular, we will reformulate the approach used in [27] in terms of bosonic fields, rather 
than in terms of Majorana fermionic fields. Our formalism allows for a direct derivation of the subleading, finite 
temperature/frequency corrections to the fixed-point dynamics. 

In bosonic coordinates, the Kondo interaction Hamiltonian is given by: 

H]? K = J^Sj : e^ P (0) cos( * sf(0)) : +S - : e -i*. P (0) :: cos (* sf ( )) : } 

= JS d -[E A (0) + E B (0)] , (14) 
where the spin densities ^a/b( x ) are given by: 

& A/B (x) = ~ [* sp + / - **](*)] ; 2% B {x) = ±= : e ±*[*.p+/-*. f ]W : . (15) 
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In Appendix A we prove that both T,a(x) and T,b{x) are SU{2) spin-1/2 operators, and show that the corresponding 
spinors at a point x are realized as: 

\a,A) x =: e^p+^K*) : |bvac) , 

and: 

1(7,5), =: e'* [^p-*"]^) : |bvac) . (16) 
The doublet \a, A/B) x provides a spinor representation of the SU{2) group generated by 

^a/b = J dyf, A / B {y) . (17) 

Also, we obtain 

Z A \o,B) x = Z B \a,A) x = , (18) 

because 

[g A)e ±i[* sp +* sf ](x)] = [g B)e ±f[*sp-*sf](x)] = _ 

Eq.(18) states that, if at a point x the spin density associated to Y*a is 7^ 0, then, at the same point, the spin density 
associated to £_b is = 0, and vice versa. 

Such a statement is the key argument used in Ref. [27] to argue that, within such a regularization scheme, the 
finite-coupling fixed point is actually moved to an infinite-coupling point. The argument is that, since it is not 
possible to have at the same point both spin densities different from 0, it is also impossible to produce a more-than- 
1/2-spin composite at the origin to overscreen the impurity spin. Therefore, the unstable overscreened fixed point 
disappears and NFL-behavior is reached at an infinitely-strongly coupled fixed point, where the impurity spin will be 
fully screened in a localized spin singlet. Such a singlet must be formed either between Sd and £a(0), or between 
Sd and Ss(0). Therefore, at the fixed point the system can lie within either one of the two singlets |Sin, A, {£}), 
|Sin,B,{S}>, given by: 

1 

ISin, A, {-}) = {^j 4 ±{\ ft) ® I I, A, {-}) - I ft) ® I T, A, {-})} 



|Sin,5{~}) = 



V L J V2 



{\1l)®\l,B,{~})-\il)(g>\tB,{~}}} 



(19) 



where 77 is the convergence factor (see Appendix A for details) and | f|~) , | JJ.) are the two impurity states with opposite 
spin polarizations. | |, A/B, {£}), and | J., A/B, {S}) are states of conduction electrons with an f or a J. A/B particle 
at x = 0, respectively, the state of all the other delocalized particles (globally denoted by {S}) being unspecified for 
the time being. 

Let us, now, define the operator Q x : 



/2ttt] 



-i"9 s f (x) 



(20) 



It is straightforward to verify that [Q x ] 2 = 1, Vx. Therefore, its eigenvalues are ±1. 

Since [Q , Hjf K ] — 0, where Q = Q x =a one may build simultaneous eigenstates of the two operators, given by: 



|Sin,+,{S}) = ^=[|Sin,A,{S})+z|Sin,B,{S})] 



|Sin,-,{S}) = i=[|Sin,A,{S})- J |Sin, J B,{S})] . (21) 

Notice that these states are not eigenstates of the spin-flavor. 

The construction will be repeated in the next Section, where we switch back to the fermion representation. 
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IV. FIXED POINT FERMIONIC S MATRIX 



In this Section we reformulate the results of Section III in fcrmionic coordinates. We will eventually get to the 
formula for the appropriate scattering potential in the unitary limit. Finally we will take the continuum limit of our 
result, and derive the single-particle S — matrix for scattering off the impurity. 

In bosonic coordinates, the two different representations for the lattice fermionic fields c aa (x) 7 which we will refer 
to as cJa 1 (x), are given by: 



JJ/ii 



(x) 



/2ttt7 



: e 2 



[*ch(x)+o-* sp (x)+a* s f(x)+/-a<7* s f(x)] . 



(22) 



In terms of cJa 1 (x), one may write down the relevant operators acting on one-particle states: the identity operator 
at point x 



^[ C ^(x)cL(x)+4^(x)c^(x)] 



(23) 



and: 



Q x =iJ2(™Mi(x)c I J a (x) - cyj{x)cUx)] 



(24) 



At a point x, Q x swaps representations / and II with each other. Moreover, since in bosonic coordinates, at x = 0, 
[Qo, H]f K ] = 0, we require the same thing to hold in fermionic coordinates. 

Following Nozieres' approach, we construct an effective fixed point Hamiltonian by introducing an infinite-strength 
repulsive potential scattering at the origin, and by making it commute with Qo- Therefore, it is given by: 



V fp = To { Um [A ]T 4" (0)O0)]P 



(25) 



where Vq — ^[lo + 2oL an d A is the strength of the interaction. 

We obtain the unitary Hamiltonian by adding the lattice kinetic energy term to Vf p . By taking the continuum limit 
of the corresponding operator, one gets the "Nozieres like" fixed point Hamiltonian. At finite A, this is given by: 



J rrr. ^ 



1 

^Tx I o i 



+ X6(x)- 



i(aa) 



-i(aa) 
1 



(26) 



This Hamiltonian envisages a scattering process for the even component of the field, which is represented pictorially 
by the boxed field in Fig.(l). 




r=lxl 



\|/(r) 

-k,l 



FIG. 1. Pictorial sketch of the even-odd decomposition of the scattering process in the (x,t) plane [35]. The heavy dot at 
the center is the scattering center located at x = 0. Ingoing and outgoing single particle wavefunctions are shown. One of the 
chiral components in the scattering has been boxed. 
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At x = the projected scattering potential exchanges the two representations. Accordingly, Hf CK is diagonal 
in the representation index everywhere, except at the origin. In the following, we use H 2 X CK to derive the fixed 
point one-particle Green's functions. From the Green's functions we derive the fixed point S matrix, which we will 
eventually use to compute the fixed point conductance. 

Let us write down the equations of motion for imaginary time one-particle Green's functions, derived from the chiral 
Hamiltonian of Eq.(26). For simplicity, here we just consider a =| and a = 1. The Green's functions are defined as 
(a, b = 1,11): 



G^(x,t;x ,t) =9(t-t')Tt 



P -0H 

-—^{x,t)^\x',t) 



-9{t -r)Tr 



(27) 



where [3 — 1/ksT and Z is the partition function. 

The coupled equations of motion for G 1 ^ and G 1 ^' 1 read: 



d_ 



IVf 



dx) 



Gl'i(x, t;x ,t ) = S(t — t )S(x — x ) — XS(x)[Gi'^ (x,t;x ,t) + iGt 1 ^ 1 (x, r; x , r )] , 



ill J, 



and: 



d d\ 
- +Wf d-x) 



dr 



G I ^ I 1 ' I (x, t;x , t ) = — A<5(;r)[Gj{' J (x, r; x , r ) — iG 1 ^ (x, t;x ,t )] . 



iii, i 



(28) 



(29) 



Eqs. (28,29) have to be supplemented with the set of equations of motion for G|{' 7/ and for Gjj^ 1 , with the replace- 
ment i — ► — i. 

The detailed solution of Eqs. (28,29) is provided in Appendix B. The S'-matrix will come out to be diagonal with 
respect to any index but the representation index, according to the definition: 



G a a £(iuj m ;x >0>x) = ^S^(iw m )G c a t(iuj m ;0> x > x) 



(30) 



where oj m are the Matsubara frequencies. 
From the results of Appendix B, we find: 



Sl'i(iw m ) = l + 2m 



A 



1 + 2\T(iu m ) 
A 



V - 27r i + 2A^ m ) 



b n ' (iu! m ) = S 1 -[ {iui m ) , b n (iu> m ) = -b n [iu) m ) 
Going back to real time, we get: 



T{iuj m — > uj + if]) = —iir 1 In 

' 7T 



D — to 
D + u) 



where D is a band cutoff energy. 
Finally : 



1 — 2Aln 


D-ui 
D+ui 




1 - 27riA — 2A In 


D-ui 
D+u 



(31) 
(32) 



S" J (u;) = -S^"(u) = 



2ttA 



1 - 27r«A - 2Aln 



(33) 
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This S— matrix is perfectly unitary in the representation space, the unitary limit is achieved with A — > oo. In this 
limit, the on-shcll S matrix in the (/, II) space, S aa {uj = 0), is given by: 

According to the definition of phase shift 5 for elastic scattering, S — e 2lS , we obtain the corresponding phase shift 
in the various channels, in the unitarity limit, given by (5| 1 // (w = 0) = — f , 1 (u = 0) = \ . In the one channel 
spin— 1/2 Kondo effect, the phase shift is 5 a = §• 

The appearance of matrix elements that are off-diagonal in the representation index avoids the unitarity paradox, 
at the price of alleging scattering processes swapping the representation, corresponding to changing the spin-flavor 
quantum number by ±2. 

The derivation above provides also the boundary conditions for the bosonic fields SS?x (A = ch, sp, fl, sf) at the 
origin. As x crosses the location of the scattering center (x — 0), none of the field changes, but \I/ S f, according to: 

* rf (a;) -» -* sf (a;)+7r . (35) 

In the following Sections, we abandon the /, //-representation, and define the physical conduction states in fcrmionic 
coordinates. To do so, we must derive the irrelevant operators providing the leading corrections to the T=0 limit of 
the scattering dynamics. 



V. THE FIRST IRRELEVANT CORRECTION TO THE FIXED-POINT HAMILTONIAN, PHYSICAL 
STATES AND THE UNITARY LIMIT OF THE CONDUCTANCE 

In the previous Sections, when bosonizing the fermionic fields, we stressed the representation redundancy associated 
to the spin-flavor quantum number. Such a redundancy does not affect the physical validity of the bosonization 
procedure, as both representations possess the same observable quantum numbers. However, it must be taken care 
of, somehow, when characterizing the fermionic Fock space for delocalized particles. To recognize what the physical 
fermionic conduction states are, we have to discuss the first irrelevant correction to the fixed point Hamiltonian H 2CK 
in Eq.(26). In order to do so, we will first go back to bosonic coordinates, and then derive the first irrelevant correction 
to the bosonic fixed point Hamiltonian Hff K of eq.(4), by building the perturbation theory about the fixed point, 
that is, by assuming the boundary conditions in Eq.(35), and by applying the Schrieffcr- Wolff procedure to compute 
the various operators. 

Irrelevant correction arise when the kinetic energy term is added to H]f K , and high energy states with a localized 
triplet at the impurity are alleged to take part to the scattering process, as virtual states (each singlet state has an 
energy Es = — 9J/4, and each triplet one has energy Et — 7 J/4, see Appendix C for details). 

To start the derivation, let us notice that, as we show in Eq.(14), the Kondo interaction Hamiltonian only contains 
spin and spin-flavor bosonic fields. Therefore, we may factorize out both charge and flavor fields, and write the kinetic 
energy in the "reduced" bosonized form: 

£ (IN 

A=sp,sf 

This factorization might be thought of as an artifact of the long-wavelength expansion, and it may be possible that, 
during the rcnormalization procedure, some terms arise, coupling charge and flavor to the remaning degrees of freedom. 
However, we will assume that such terms are irrelevant anyways, at low enough temperature. 

Using the Schrieffer- Wolff procedure, we take as the lowest-energy subspacc the one spanned by the singlets 
| Sin, A, {£}), | Sin, B, {S}), and we construct an effective Hamiltonian as a perturbative expansion in t 2 /J. Defining 
the projector onto the lowest-energy subspace as : 

P = ]T |Sin, «,{-}) (Sin, , (37) 

u=A,B 

the effective Hamiltonian, up to terms 0(t 3 /J 2 ), is: 

#Red[l - Po]ZfRed 



(36) 



i/ E ff « P (ffRcd + H K )P Q + P 



Eg — Ej 
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1 ^ 2r 



Eg — Ej* 



#Red[l - Po]#Rcd[l - Po]-ffRod 



Po ■ (38) 



In this Section we focus onto the first term at the r.h.s. of Eq.(38), Po(//R C d + Hk)Pq- In Appendix C, we derive 
the action of H^cd hi the discrete lattice model, where H^cd is substituted by the corresponding lattice operator, 
Hr,Red (see Appendix C for details). When computing matrix elements of Ht, Red between different impurity states, 
we include only the impurity neighboring sites (which is equivalent in spirit to Wilson's NRG approach), and use the 
symbol hr, when referring to the corresponding operator. We show that, once projected on the space of the singlets, 
the corresponding hopping term takes the form: 

P H Red P ^ ~tQo(Q a + Q- a )Qo ■ (39) 

When written in terms of fermionic fields, the operator in Eq.(39) contains contributions that are off-diagonal in 
the /, II representation. However, in the continuum limit a — > 0, such terms just add up to the scattering potential 
in Eq.(26), so that they can accounted for by substituting A — > A — t. This has no consequence on the developements 
of Section IV, because A — > oo at the end, but shows that states / and II have to be properly mixed by hopping at 
any distance from the impurity. 

We recognize that the corresponding physical requirement is that Q x commutes with the Hamiltonian not just at 
the origin, but at any point x. Hence, physical states can be constructed by using the projection operators given by: 

v± = ^H(i x ±Q x ) (40) 

* X 

and by requiring that pyhsical one-particle states are unaffected under, for instance, application of V+. These states 
are given by: 

|phys+, aa) x = -±= [cli(x) - i^a)^ 1 {x)] |0) (41) 

(The other set of one-particle physical states, |phys— , <ra) x , is obtained by changing i — > — i). 

Therefore, the transformation matrix U, that maps /, //-representation onto the basis of the physical states, is: 

\ l/s/2 i/V2 J ■ 
U transforms the S matrix of Eq.(34) as follows: 

uts - u =("o 1 2) • < 42 ) 

As we have projected onto just one of the two possible sets of physical even symmetry states, we find an S matrix 
for the even states which has just one non vanishing element (that is -1, implying that also the even wavefunction 
has a node at the origin, as it must be in the unitary limit). The other possible physical states decouple from the 
scattering dynamics. Had we chosen the other set of states, the transformed matrix would be, instead, with only the 
bottom diagonal element different from 0. 

At this point we have to recall that this result refers to the even-parity wave. In fact, we get full transmission for 
both parities I = e,o, because the ^-matrix is diagonal in the basis of physical states, and the matrix elements are 
given by S^'' = —1, for both parities I = (e, 0). The transmission across the impurity is given by the trace [36]: 



T = Tr 



(+) 



Tr (+) {1( . (43) 



Here Tr^ + j means tracing over all physical degrees of freedom a, a. Although the system posseses 4 degrees of freedom, 
they are shared by the two physical zb-ground states. In fact, there is a Z2 symmetry breaking in this system. Hence 
7> (+) {l} = 2. 

By applying Landauer formula to Eq.(43), we obtain the correct conductivity in the unitary limit: 
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2e^ 
h 



(44) 



This corresponds to the halving of the zero-point entropy of the ground state, which turns to be l/21n(2) for each 
ground state. 

The fractionalization of the degrees of freedom can be understood also within the (I, II) representation for the 
S— matrix given in Eq.(34). One can think of closing the line onto itself in a symmetrical ring geometry, with two 
equal impurities at opposite sites. In this representation, in the upper branch we have propagation of type-/ forward 
scattering into type II across one impurity, while in the lower branch the scattering is back, 77 — ► /, across the other 
impurity. Each process gives rise to a phase shift at uo — 0, given by S e — ±7r/4. The conductance across each one 

~ 2 o 2 

impurity is G = 2 x 2 x ^- sin S e = 2^-, where the first factor comes from the spin degrees of freedom, while the 
second one comes from the number of channels. Thus, the unitary limit is again obtained. 



VI. LEADING FINITE-T CORRECTIONS TO THE UNITARY CONDUCTANCE 

In this Section we will explicitly write down finite-T corrections to the fixed point conductance, coming from higher- 
order corrections to the fixed point Hamiltonian, (that is, from operators arising in the perturbative expansion of 
Eq.(38)). In Appendix C we derive term by term the contributions to Eq.(38), ut ro third order. The matrix elements 
of the various operators in the basis of the singlets |Sin, u, {£}), (u = A, B) form a 2x2 operator matrix, acting on 
the two low-energy singlets. Eventually, we will restrict them to the subset of physical states defined in the previous 
Section. 

In particular, in this Section we will derive finite-frequency contributions to the fcrmionic 5-matrix. These will 
come out to be oc y/uj, which shows the NFL- nature of the corresponding ground state [11]. Vertex corrections provide 
higher-order contributions, which we will not consider here. 

Following the derivation of Appendix C, we see that, to 0(t 2 /J), we obtain the following matrix elements: 

Of 2 f 2 

M 2 AB — M BA = ; M\ A = M% B = ^-jj « -— . (45) 

These terms provide just an over-all trivial shift of each energy eigenvalue by a constant amount 0{t 2 ). 
Non trivial effects, instead, arise from the third-order corrections. 
From the calculations reported in Appendix C, we see that: 

M\ A = M BB = ; M\ B = sin [* sf ( )]^|^ (46) 

plus terms 0{t 2 / J) (and higher) that renormalize the ones considered before. 

The third-order correction in Eq.(46) is not affected by the physicality constraint because of the operatorial relation: 
\Y[ X 2xj sin[^ s f (x)]] — 0. Since M 3 changes by ±1 the spin-flavor, the corresponding diagonal (in /, //) contributions 
to the fixed-point Green's functions will be zero. Instead, it gives an 0(t 3 /J 2 ) off-diagonal correction to G 1 ^ 1 . In 
the interaction representation it reads: 

6G lZ' a '( X > T '> X '> T '^ = r^Tvje-^T^^Ja;,^^^)^ 7 ^'^')]}^^ . (47) 

J 

Computing the trace of Eq.(47) requires using the bosonic representation for the operators (fi 1 ' 11 , provided in Eq.(22). 
In bosonic coordinates, Wi nt (ri) is defined as: 

Wmt(ri) - ^ eTlH!p i- sin[* rf (0)] = ^W^^ ■ ( 48 ) 

Two independent periods appear in the Green's functions: the length L and the inverse temperature (3. A rigorous 
calculation of the correlators in Eq.(47) would require the introduction of Jacobi's elliptic ^-functions. However, here 
we will attempt to compute just the finite- frequency corrections, at T = 0. Therefore, we shall approximate Eq.(47) 
as: 
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t 3 

irJ 2 



|* ch (£c,r) 



:>(: 



-ai*n(i,T) .. aj*fi(z ,r ) 



(:e 



-<75* sp (x,r) 



9*2 

9a; 



■(0,n) : e CT 2* sP 



}(:e 



+ffajt at (i.i) .. 



:: sin[* sf (0,n)] :: e 



-cra^^'sf (a; ,r ) 



) , 



(49) 



where (. . .) denotes ground state average. As shown in Appendix D, for the case {aa} — {f 1} and in the limit of 
T = and large L, Eq.(49) provides the result: 



SG 1 ^ 11 (iui m ;x,x ) = 



e / 



J 2 v] V nv f 



(50) 



with x — x' > 0. In Eq.(50), we have performed the analytic continuation to real frequencies for the retarded Green's 
function: Lo m — > — iu>. 

Introducing in Eq.(50) the Kondo temperature as the relevant physical energy scale according to the substitution: 
2^3/2 VfJ 2 — * we obtain a leading finite- frequency correction that goes as w5, in agreement with the results 

obtained in [11]. In particular, we have found a finite-w correction to the S matrix given by: 







-i{aa)J^ 







By projectig the result of Eq.(51) on the basis of physical states, we get: 

Ut[S ffQ +<5S CTQ ]U = 



(51) 



-1 + ,/H. 











(52) 



As a side remark, let us consider the case of a QD hybridized to metal contacts by a tunneling potential V. If we 
use the relations found within the Anderson model [38], we find the width of the Kondo resonance to be given by: 
r = 4kB J 1 ' { = 7ri^(0) I 1 2 , where V is the tunneling strength and v(0) — 2tt/vj is the density of states at the Fermi 

level. This implies, in our case, that |V| = ^-J2-kv^. 



By using the result in Eq.(52), we compute the finite-energy transmission: 

2 ' 



T(w,T = 0) =Tr {+) 



2 1 



to 



(53) 



Finally, in order to obtain from Eq.(53) finite-temperature dependence of the conductance, we have to recall that, 
at finite temperatures, there are two contributions to the conductance: one arising from the smearing of the Fermi 
surface, the other from the inelastic processes. In fact, for the large U Anderson model, the explicit dependence of 
the transport time on temperature at the Fermi energy has to be taken into account separately [38]. This gives: 



em = e jJ 



dw 



duo 



T(w,T) 



2e 2 r°° 



+ e^) 2 



1-1 




2 " 


2 


T K 





+ A(/i,T) 



2e^ 
h 




(54) 



where f(u>) is Fermi distribution, and A is the total contribution coming from inelastic processes, which we neglect 
here, as it is assumed to provide corrections that are higher order than \[T . Eq.(54) contains the ultimate result of 
our derivation: the calculation of the fixed-point contribution to the conductance, together with the leading finite-T 
correction, and the elucidation of the connection between this correction and the various scattering processes that 
take place at the fixed point. As we already mentioned, our formalism allows also for calculating inelastic term arising 
from vertex corrections, but we will not consider them here. 
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VII. CONCLUSIONS 



In this paper, we use Landauer formula to derive the conductance at the two channel spin 1/2 overscreened Kondo 
fixed point, together with leading finite-temperature corrections. We perform our calculations within the framework 
of our simple Hamiltonian theory, in which we derive a suitable fermionic representation of the S— matrix, by using 
bosonization as an intermediate step. In order for us to achieve such apparently simple results, we had to go through 
various mathematical approaches, which are used in the literature to analyze different aspects of the overscreened 
Kondo problem [11,27,34]. For instance, we had to complement the regularization scheme introduced in Ref. [27] 
with the bosononizaton technique widely used in [11], and with the careful discussion about the role of the spin-flavor 
quantum number in Ref. [34]. 

Our research is motivated by the renewed interest in the Kondo model, recently arisen in connection to conductance 
experiments across QD's. In general, using simple models for correlated electros, like the Anderson model, allows 
for grasping the physics involved in tunneling experiments across confined areas between two contacts, as well as 
across Coulomb blockaded systems (like the QD device we have in mind) . Recently, efforts have been made to achieve 
an exact description of transport also in the nonequilibrium case, starting from the intcgrability of the two-lead 
Anderson model [37], which confirms previous numerical RG results [23]. Despite the exactness of these results, they 
are unsuitable to our case, as they refer to the Fermi liquid 1CK fixed point. 

Recently, various groups have been predicting that 2CK could be realized in QD systems. In particular, we have 
proposed that orbital 2CK effect may arise in a vertical structure with cylindrical symmetry around an applied 
magnetic field [16]. 

In our proposal, the dynamical degrees of freedom involved in the scattering across the interacting Dot are provided 
by an appropriate combination of the transverse angular momentum of lead electrons, m, and of their spin, a. The 
single particle wavefunction of the delocalized electrons traveling along the z direction has an orbital part factorized 
in the cross-sectional plane. In a linearized band picture and close to the Fermi energy cf, their wavefunction of 
energy e and parity I w.r.to z = is : 

where 5 = e — hvp\q\, and \a is the spin wavefunction. 

Another proposal is based on the use of an additional lateral QD to tune the exchange coupling to the channels 
to the isotropic point [17]. In such a geometry, however, even a small anisotropy in the coupling to the channels is 
enough to drive, at low enough T, the system from the 2CK fixed point to a 1CK channel fixed point, with strong 
coupling only in the dominant channel [39]. This problem does not arise in [16], because of the assumed cylindrical 
symmetry. This symmetry does not allow for any off-diagonal coupling mixing the two channels, because it enforces 
the angular momentum selection rule in the cross plane. Of course, it is very demanding to produce such a strictly 
cylindrical system experimentally. 

These proposals have triggered a renewed interest in 2CK transport. On the theoretical side, reconsideration of 
the matter is relevant, in view of the particular kind of devices involved in the proposed experiments. Indeed, the 
scattering approach used in Ref. [11], in conjunction with CFT techniques, is better suited for the single impurity 
s— wave scattering in a three-dimensional medium than for a two-lead device. In [18], the authors identify the two- 
channel Kondo fixed point as a quantum critical point between two Fermi liquid phases, and derive its dependance 
on external parameters such as temperature, magnetic field and voltage bias accordingly. 

For the purpose of understanding the physics of the transport across mesoscopic devices, we believe our approach 
to be a straight connection between the physics of the overscreened Kondo problem at the fixed point and its physical 
consequences. In a clear and easy-to-follow framework that extends Nozieres' hypothesis, it sheds light on the physical 
processes that happen at the impurity as T — > 0, and relates them to the macroscopically detectable non- Fermi liquid 
behavior in the T-dependence of the conductance. The GS at the NFL fixed point is found to be degenerate. Because 
the physical system only involves one of these GS's and its corresponding excitations, the contribution of each degree 
of freedom to the conductance is halved. As discussed in Section IV, this leads to the same unitary limit of the 
conductance as in the one channel case. Finally, our approach might also provide an alternative route to investigate 
the effect of anisotropy on the conductance. 

We acknowledge fruitful discussions with M. Fabrizio, P. Wiegmann and I. Aleiner. Work partially supported by 
TMR Project, contract FMRX-CT98-0180. 
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APPENDIX A: BASIC BOSONIZATION STEPS 



In this Appendix we will review some basic bosonization steps, which will allow us to write down fermionic fields 
in terms of bosonic operators. 

In order to write fermionic fields in bosonic coordinates, let us introduce a massless scalar bosonic field ^(x), given 
by: 



n#0 



n 



where rj = + is a regularizator. 

The algebra of the bosonic modes is: 

[<l,p]=i ; [a n ,a m ] = n5 n+m fi ■ 
Therefore, we may split ^(x) into a creation and an annihilation part, ^(x) — ^+(x) + ty_(x), with: 

^+{x) = q-iy e l x L n ■ y_f x \ = + , 

^— ; n L 

n— 1 n— 1 

and: 

[*-(!/), *+(a;)] = - ln(e 2 * a,/ - e 33 *?* (*+*")) . 
The ground state of the Fock space spanned by the bosonic modes, |bvac), is defined by: 

p|bvac) = a„|bvac) = (n > 0) . 

A fermionic field may be defined as: 

c\x) = :e 4 *^ : 

where the columns : : denote normal ordering with respect to |bvac). 
Indeed, by using the general identity 

e A e B = e B e A e^ , 

which holds if [A, B] is a number, one gets: 



. e iat>(x) .. e -i/3*(y) ._ gi[a*+ (x)-/3*+(j/)] gi[a*_ (x)-/3*_ (j/)] 



As a, /3 = ±1 and x ^ y, we get: 

c(ai)c(y) + c(y)c(x) = c(x)c'(y) + c'(y)c(x) = 



(Al) 



(A2) 

(A3) 

(A4) 
(A5) 
(A6) 



(A7) 



(A8) 



On the other hand, by properly regularizing the anticommutator we get, as x — > y: 

while as x — > y the anticommutator between c(x) and c(y) remains equal to (as well as the anticommutator 
between c^(x) and c^(y)), one gets: 



c{x)c\y) + c'(y)c(x) = 5 



x-y 



(A9) 



The density operator at a point x may be expressed in bosonic coordinates by means of the "point-splitting 
regularization" as follows: 



p(x) = lim : {: e**^ :: e"^ :} := _1^£M 



(A10) 



The free kinetic Hamiltonian may also be written in a bosonized form as follows: 
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H T = -iv f J dx : e»+W e '*-Wl[ e - , »+We- it -W] := J dx 



dx 



(All) 



When fermions carry several quantum numbers (spin, flavor), bosonizing requires first of all introducing many 
bosonic fields Also, in order to make fermions with different quantum numbers anticommute, one has to introduce 
a "Klein factor" r] aa in front of each bosonized field. In general, one chooses r\ aa to be a real Majorana Fermion, that 
is, {r\ a ,a) 1 = 1 [19]. As a consequence, fermionic operators defined as: 

Caa{x) = rj aa : e -5[*ch(x)+o-* sp (x)+a*fi(x)+aCT* sf (x)] . . c f ^ _ ^ . g f [* c h(x)+CT* sp (x)+a* fl (x)+aCT* s f (x)] . (A12) 

anticommute with each other for different quantum numbers, as they must do. 
The following commutation relation holds: 

rJ_d*(x) j:eja * (B) , = , :e<a *( B ) . _ (A13) 

2tt dx 

By using only one bosonic field &(x) it is possible to build an SU{2) spin-current j(x). The components of the vector 
current density are given by: 

W-^^iM-iSd). (A14) 
Indeed, according to Eqs.(A7), in the case a = (3 = \[2, and to Eq.(A13), we obtain: 

[j+{x),j-(y)\ = i^\x -y) + 5(x- y)f{y) (A15) 

and 

[f(x),j ± (y)]=±5(x-y)j ± (y) . (A16) 

Eqs.(A15,A16) provide us with the usual SU(2) affine algebra obeyed by the spin-density operator. 
This, in particular, proves that the operators Ha(x) and Eb(i) we defined in Section III are SU(2) spin current 
operators. 

The corresponding spinors at a point x may be created by acting |bvac) with : e ± 5[* B p+* rf K a; ) : anc | : e ±f[*sp-*sf](x) 
respectively. Let us define: 

\a;A) x = [: e^p+'MM :]|bvac) ; \a; B) x = [: e a ^"- 9 "^ :]|bvac) . (A17) 

The following commutation relations hold (either the upper or the lower signs hold): 

\Ei(x), : e ± 5 [*.p+*-](») : ] = [£±(x), : e ±il*.P-*-](w) : ] = (A18) 



and: 



[£±(x), : e =Fi[*.P+*.f](») : ] = S (x - y) : e ±i[*.p+*-K») : ; 



[S|(x), : e T 5[*sp-*sf](!/) .] = S (x - y) : e ±§[* sP -* sf ]fe) . . ( A 19) 



Finally: 



[£^(a;), : e ± *[*'«' + *' f ]^ :] = ±-S(x - y) : e ±§[* B p+* S fl(i/) . . 



2 



[£fj(x),: e ± 5[*=p-*= f Kf) :] = ±i$(a; - y) : e^*^"*^) : . (A20) 

The set of equations listed above shows that the doublet \a, A) provides a spinor representation of the SU(2) group 
generated by Sa, and that the doublet \a,B) provides a spinor representation of the SU{2) group generated by 
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To conclude this Appendix, let us now prove that fermionic fields belonging to two different representations anti- 
commute with each other. 

Let us start from the fields in the two representations expressed in bosonic coordinates: 



4>l a i x ) = V™ ■ e 



-^[*ch(a;)+<T*3p(a;)+a*fi(a;)+a<7* B f(a;)] . 



and: 



^N aa - i[* ch (x)+<7* sp (x)+a* fl (x)-aCT* sf (x)] . 



where the Klein factors r\ aa are given by real Majorana variables. 
Clearly, fields within the same representation anticommute: 

{tf a (x),<l>yAy)} = 5 t „ T .5 aa .5(x-y) . 



(A21) 



(A22) 



(A23) 



Let us, now, consider fields from the two different representations. Let us start with two annihilation field operators 
(in the case {aa} ^ {<r a }) : 



}} ■ exp 



--[*ch(a;) + *ch(2/) + CT*sp(s) + o * sp (j/)- 



a* fl (x) + a * fl (y) - aa^ s{ (x) - aa * s f(j/)] 



:= 



By following the same procedure, we obtain the result: 

{<i> I J a (x),<t> I Jl(y)} = 

and the anticommutation relations: 

{4>l a (x),<f> I Jl,{y)}= Vaa r, a . a ,{eW 1+ ™ 



(A24) 



(A25) 



[6-^ _ e ^]"i[l+- +aa — aa ] ]} . cxp[ _ ^ ch{x) _ ^(y) + ^ _ CT '* 8p (j,) 



and: 



a* fl (x) - a * fl (j/) - acr^ s f (x) + a a % f (y)]] := , 



W I J a (x),<j> I J a ,(y)} = o , 

that complete the proof of the identities used throughout the paper. 



(A26) 



(A27) 



APPENDIX B: CALCULATION OF THE FIXED POINT GREEN'S FUNCTIONS 



In this Appendix we will show in detail how to calculate the fixed point Green's functions with the method of the 
equations of motion. Within our framework, this will come out to be a straightforward application of one-dimensional 
scattering theory. 

The equations of motion for and G|{' 7 are reported in Section IV and are given by: 

(-^ - iv f^~\ G\i(x,t; x ,t ) = 8(t - t )S(x - x ) - XS(x)[G^ (x, t;x , r ) + iG 1 ^' 1 (x, t;x , t )] (Bl) 
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and: 



d_ 

97 



ivf-^-^j G 1 ^' 1 (x, t;x , t ) = — A<5(x)[G|{' / (x, t; x ,t ) — iG^(x, t;x ,t )] . 
To solve the set of Eqs.(Bl,B2), let us introduce the Green's functions in the mixed representation: 
Gf 1 B (x,T;x',T) = ^J2[ dpe-^e^G^(iu miP -x) . 

ioj m 

In the mixed representation, the equations of motion become: 



(B2) 



(B3) 



{-iw m + v f p)G^(iuj m ,p;x') = e tpx - X J dq [G 1 ^ (iuj m ,q; x' ) + iG" J (iuj m , q; x )} (B4) 



and: 



(-iuj m + v f p)G I l I { 1 (iu; m , p; x' ) = -A J dq [G" J \iu m ,q:x ') - iG\{ (iuj m , q; x ')] 
The solution to the set of Eqs.(B4 , B5) is given by: 



(B5) 



G^\ (iLU m ,p;x ) = 



-tpx 



lU) m + v f p -iu m + V fp 



X 



1 + 2X'T{iuj m ) 



/ p -iqx \ 
dq — (B6) 



and by: 



G n ' (iuj m ,p;x ) = — 



iA 



-iw m + vfp [1 + 2A' ' T(i<jj m ) 
where A = A/v/ and T{iuj m ) is defined as: 

F{iuj m ) = v f { dp — — 

J -iu m + v f p 



J dqv } (- 



iw m + v f q 



In 



D + iu> ri 



-D + iu r 



(B7) 



(B8) 



In Eq.(B8), D is a high-energy band cutoff. 

In order to derive the S'-matrix elements, in Section IV we use the real space Green's functions for x > and x < 
in the limit A — > oo. In the case uj m > they are 



G^l{i(jj m ;x,x ) = J dpe vpx G^l{iijj m ,p]x ) = 



1 - 2m 



1 



2T{lUJ m ) 



and: 



G"' 7 (iw m ;x,a;') = f dp e^Gl 1 / (iu m ,p; x) = — 6{w m ) 

J Vf 



2m 



2T{iu m ) 



The "noninteracting" Green's functions Gi°/' a ' b (iuj m ; x, x ') for x — x > are given by: 



^(0)50,6/ . 

G\{ (iu m ;x,x 



dp 



e i P (x-x ) 2m 



Vfp - iw m Vf 
From Eq.(Bll), we see that Eqs.(B9,B10) take the form: 



6(w m )e "7 {x x) 5 a > b 



G^(iu; m ;x,x') = Gf}' JJ {iuj m ]x,x) 



and: 



G^(iuj m ;x,x) = G^'' 1 ' 1 ' (iu m ;x,x) 



(0);/,// 



1 - 2ni 



2m- 



X 



1 + 2X'J 7 {iLu m ) 



iX 



'l + 2X'T(iu m ) 

From Eqs.(B12,B13), we derive the S matrix elements in Section IV. 



(B9) 



(BIO) 



(Bll) 



(B12) 



(B13) 
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APPENDIX C: ACTION OF H T ON THE VARIOUS STATES 



In this Appendix, we show in detail how hr acts on the states arising from hybridization of the impurity spin with 
the spin of itinerant lead electrons. This is the mathematical support to the derivation of Section V, where we use 
the results of this Appendix to derive leading corrections to the fixed-point Hamiltonian. 

The triplet states are: 

A, 1, {H}> = | T, A, {H}> ® | -fr) ; |Tri,S,l,{S}} = | T,B,{S})®|^) ; 
|Tri, A, -1, {S}) = \ i,A, {S}) ® | JJ-) ; |Tri, B, —1, {£}) = | [,B, {S}) ® | JJ.) ; 

|Tri,A0,{S}) = T,A{5»®U> + | l,A,{E})®\ ft)] ; 

|Tri,B,0,{S}) = i=[|T,S,{H})®|^ + ||,B,{S})®|^] . (CI) 
Let us introduce the lattice operators bx given by: 



b x (x) = ^/^: e-**™ : ; X = sp,sf . (C2) 

To start the derivation, let us rewrite H Rcd in a lattice form, by using the operators b sp , 6fl. We get the reduced 
lattice kinetic energy operator H T ^ ed , given by: 

^T,Red= ^ {6^(a;)6x(a:)-t[6^(a;)(6x(a: + a) + 6 x (a:-a)) + (&^(a: + a)+6^(x-a))6 x (a:)]} . (C3) 

a;,X— sp,sf 

From the operator in Eq.(C3), we may single out the term that acts on the Fock spaces of the states formed by 
hybridization between the spin of the localized impurity at the origin, and the spin of conduction electrons at x = ±a. 
Such a term is given by 

hr = -t ]T [b x (0)(bx(a) + b x (-a)) + (b x (a) + b x (-a))bx(0)] . 

X— sp,sf 

The action of hr on the singlets is given by: 

HSin, A) = -j={[b sp (a) + M-a)]|Tri, B, 1) - [6|» + &| p (-a)]|Tri, B, -1)} 



-^[6 sf (a) +M-a) + bl(a) + ?4(-a)]|Sm, B) + ^[6 sf (a) + 6 sf (-a) b^ (a) - 6^ (-a)]|Tri, B, 0) (C4) 



and by: 



ftrlSin, fl) = -±={[b sp (a) + 6 sp (-a)]|Tri, A, 1) - [&t p ( a ) + &t p (-a)]|Tri, A,—l)} 



— Ma) + b sf (-a) + bl(a) + (-a)]|Sin, A) - 1 [6 sf (a) + 6 sf (-a) - b^a) - ^ f (-a)]|Tri, A, 0) . (C5) 
When acting on the triplet states, instead, hr provides the following results: 

^ITri,^, 1> = -t[6+ f (a) + 6^(-a)]|Tri,,B,l) + [ & | p(a) + 6 t p( _ a)] | Sin5 ^ _ [ fo | p(a) + 6 t p( _ a)] ]Tri! 0> (C6 ) 

hr\Tn,B,l) = -t[b sl (a) + b sf (-a)]\Tri,A,l) + -^=[&| p (a) + &| p (-a)]|Sin, A) - -^=[6+ p (a) + fo+ p (-a)]|Tri, A, 0) (C7) 
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hr\Tri,A,-l) = -t[b s{ {a) + b s{ (-a)]|Tri, B, -1) - -^=[b sp (a) + 6 sp (-a)]|Sin, B) - -^=[6 sp (a) + fe sp (-a)]|Tri, S, 0) 



(C8) 



hr\Tn, B, -1) = -t[6^(a) + ^ f (-a)]|Tri, A, -1) - ^=M«) + & sp (-a)]|Sin, A) - ^=M«) + & sp (-a)]|TVi, A, 0) 



h T \Tn, A, 0) = -^={[b sp (a) + 6 sp (-a)]|Tri, B, 1) + [&|» + 6 s t p (-a)]|Tri, B, -1)} 



(C9) 



-gUM") + + &k(a) + 6t f (-a)]}| 1 &i,B,0) - [6 sf (a) + 6 rf (-o) - 6 s f f (a) - ^ f (- a )]|Sin, B» 



(CIO) 



/i T |Tri, B, 0) = -_{[6 sp ( a ) + 6 sp (- a )]|Tri, A, 1) + [&J p ( a ) + 6| p (-a)]|Tri, A, -1}} 
-{[6 sf (a) + b si (-a) + 6 s f f (a) + 6^ (-a)]}|Tri, A, 0) + [6 sf (a) + 6 sf (-a) - to - ^ (-a)]|Sin, A)} . (Cll) 



The equations listed above are all we need, in order to work out sublcading corrections to the fixed point Hamilto- 
nian. 

We apply the Schricffcr- Wolff transformation introduced in Eq.(38), by including just the hopping term between 
sites nearest neighbors of the impurity. The effective perturbative pcrturbative potential is: 



Ili, : P () |(; U,, > + 



1 



Eg — Et 



h T [l-V ]h T 



1 



Eg — E j 1 



/l T [l-P ]/l T [l-Po]'lT 



Pn . 



(C12) 



The first term at the r.h.s. of Eq.(C12), can be expressed in terms of the operators Qo, and Q± a , leading to the 
Eq.(39). 

The 0(t 2 /J) of Eq.(C12) provides the following matrix elements: 



M 2 AB - M 2 BA = ; M\ A = M 2 BB = 



2t 2 



t 2 



E-\3 2J 



(C13) 



Therefore, on |Sin, ±, {H}), second-order (in t) dynamics just yields an over-all trivial shift of each energy eigenvalue 
by a constant amount. 

Non trivial, effects, instead, arise from third-order corrections. Indeed, while we have once more: 



M AA = M% B = 



on the other hand, we obtain: 



(U) 2 M AB = J2( Si ^ A ^ E }\ h T\ E x)(Ex\h T \E x ,}(E xl \h T \Sm,B) = 

XX' 



(C14) 



3t s 



(6 sf (a) + b s! (~a) + bl f (a) + b^-a)) 



■ j[(M«) + M"*)). (bt p (a) + bt p (-a))}(b sl (a) + M"«) b^a) b^-a)) 



(C15) 
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where \Ex) is a generic high-energy triplet. 

From Eq.(C15), we see that M\ B contains a contribution proportional to the term we found to 0(t). Therefore, it 
may be accounted for by means of a slight renormalization of A , that is 0(t 3 /J 2 ) and is absolutely irrelevant. Since 
this effect is trivial, we will not consider it here. The nontrivial part of the correction operator can be derived by 
recalling the basic bosonization rules listed in Appendix A. In particular, we obtain the following expression for the 
commutator in Eq.(C15): 

[(&t p ( ffl ) + bl(-a)), (b sp (a) + b sp (-a))} = ^^g^ • (C16) 
Therefore, taking the limit a — > in the regular part of the result in Eq.(C15), we get the result: 

M\ B = M 3 = S in[* sf (0)]^g^ (C17) 
plus the irrelevant terms discussed before. 



APPENDIX D: CALCULATION OF LEADING FINITE-FREQUENCY CORRECTIONS TO THE 

5-MATRIX 

In this Appendix, we will show in detail how to calculate leading finite-frequency corrections to the S'-matrix. To 
do so, first of all, let us recall what are the relevant correlators we need. We may calculate them according to the 
basic bosonization rules of Appendix A (and to the fact that all the involved fields are chiral). We obtain: 



(T T [: e -4*< x ' T > :: e**< x ' lT '> :]) = - 



(Dl) 



e ^jr(x+iv f T) _ e ^(x'+iv f T') 

(this will provide us with the charge-charge and the flavor-flavor part of the relevant correlator) . 
Moreover, we get: 

(Tr[: e-^(^) : : e**-C T > :]) = - j—^j g- . - -L, . , • (D2) 

Finally, we obtain: 

(T T [: e ~5*sf(a;,r) .. e i* 8f (0,ri) .. e -±V 3l (x ,t) _ e f* 8f (x,T) .. e -i* s f(0,Ti) .. e §* sf (x',T') .j^ _ 



[ e — r^Ti - e 2 E i ( x + i "/' r )]5[ e — zAi - e^jrix'+ivfT)^ 



(D3) 



By putting together all the correlators in Eqs.(Dl,D2, D3), we obtain the integral in Eq.(49). To calculate the 
integral, we use as an auxiliary variable £ = exp[— ^-Vfr]. Therefore, after all the substitutions and the variable 
replacements have been made, the integral reads: 

it 3 f 2tt\^ f 1 [ e T i ( I +™/ T )- e ¥( I '+ w / T, )li 



2irvfJ 2 \LJ J q ' e 2 Tr(x+iv s r)^^_ e ^(x +iv S T')^l 

[ e ^P(x+TO/T) + g2p(x'+TO/r')] ]^{x+iv,r) + ^(x'+w/r')] _ 2 



it 3 {2w/L)i f [e^^+^fr) + e ^(x +zv f r )] [ e ^(*+ 

2-KVfJ 2 [ e ^x+iv f r) _ e 3p (l ' +i „ /T ')j3 I ^ e 3p (a . +i „ /T)e 2p (a ./ +i „ /T / ) - ^ e ^±( x + 



iv s r ) _ l)( e Tr( x +"'/ r ) - 1) 

(D4) 

(notice the extra prefactor of {2i:/L)i , which we have introduced in order to make the bosonization rules of Appendix 
A consistent with the normalization for the Green's functions introduced in Appendix B). 
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In Section VI we point out that the integral in Eq.(D4) is computed in the limit of very large L,f3. Such a limit 
yields the ultimate approximate formula: 

"' (±V ! r (D5) 



2irvfJ 2 \iJ [x-x'+ivfT]i 
When going back to (Matsubara) frequency space, we have to calculate: 



' i) / dT l V 1» ' (D6) 

1 J J-oc [X - X + IVfTl 2 



2-KVfJ 2 \i) J_ 00 T [ x -x +iv f r]i 
The integral in Eq.(D6) is calculated as follows: 

f3 / 1 \ I r°° „iumT _„'+3 



2TTVfJ 2 \lJ J_ 00 [x-X +lVfT}2 lTVfJ 2 \lJ dxJ_ 00 y/TT J_ 00 



v f) 



with x — x > 0. 

When rotating back to real times/frequencies, the result in Eq.(D7) reads: 



— ^jMe^l" ') . (D8) 

j 2 w; V 

By going back to real frequency, this provides the finite-frequency correction to the S-matrix that we have used in 
Section VI to derive the conductance. 
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